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By  utilizing  the  state  feedback  exact  linearization  approach,  a  nonlinear  robust  control  strategy  is  designed 
based  on  a  multiple-input  multiple-output  (MIMO)  dynamic  nonlinear  model  of  proton  exchange  mem¬ 
brane  fuel  cell  (PEMFC).  The  state  feedback  exact  linearization  approach  can  achieve  the  global  exact 
linearization  via  the  nonlinear  coordinate  transformation  and  the  dynamic  extension  algorithm  such 
that  robust  control  strategy  can  be  directly  utilized  to  guarantee  the  robustness  of  the  system.  The 
proposed  dynamic  nonlinear  model  is  tested  by  comparing  the  simulation  results  with  the  experimental 
data  in  Fuel  Cell  Application  Centre  in  Temasek  Polytechnic.  The  comprehensive  results  of  simulation 
manifest  that  the  dynamic  nonlinear  model  with  nonlinear  robust  control  law  has  better  transient  and 
robust  stability  when  the  vehicle  running  process  is  simulated.  The  proposed  nonlinear  robust  controller 
will  be  very  useful  to  protect  the  membrane  damage  by  keeping  the  pressure  deviations  as  small  as 
possible  during  large  disturbances  and  prolong  the  stack  life  of  PEMFC. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

As  a  renewable  energy  source,  fuel  cells  are  one  of  the  promising 
energy  technologies  for  sustainable  future  due  to  their  high  energy 
efficiency  and  environment  friendliness.  Compared  with  the  other 
types  of  fuel  cells,  a  proton  exchange  membrane  fuel  cell  (PEMFC) 
shows  promising  results  with  its  advantages  such  as  low  temper¬ 
ature,  high  power  density,  fast  response,  and  zero  emission  if  it  is 
run  with  pure  hydrogen,  and  it  is  suitable  for  use  in  portable  power 
supply,  vehicles,  and  residential  and  distributed  power  plants  [  1  -3  ]. 

PEMFC  is  a  nonlinear,  multiple-input  and  output,  and  strongly 
coupled  dynamic  system.  Its  working  process  is  accompanied  with 
liquid/vapor/gas  mixed  flow  transportation,  heat  conduction  and 
electrochemical  reaction.  The  output  current  changes  when  the 
drove  load  changes,  and  electrochemical  reaction  is  accelerated, 
simultaneously.  If  the  flow  rate  of  oxygen  is  too  low  in  cathode, 
the  output  power  of  PEMFC  system  could  be  decreased  because  of 
lacking  oxygen,  which  is  so-called  starvation.  Therefore,  in  order 
to  generate  a  reliable  and  efficient  power  response  and  prevent 
membrane  damage  as  well  as  detrimental  degradation  of  the  stack 
voltage  and  oxygen  depletion,  it  is  so  significant  to  design  an  effec¬ 
tive  control  scheme  to  achieve  optimal  air  and  hydrogen  inlet  flow 
rates. 
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At  present,  many  control  strategies  have  been  adopted  for  con¬ 
trolling  the  PEMFC  system.  Vahidia  et  al.  [4]  adopted  a  predicting 
control  method  to  design  the  controller  for  fuel  cell  vehicle  and 
it  could  satisfy  the  quickness  response  of  air  supply.  Golbert  and 
Lewind  [5]  who  made  fuel  cell  model  linearization  used  predict¬ 
ing  control  to  satisfy  the  need  of  power.  Yuan  et  al.  [6]  imposed  a 
predicting  control  method  based  on  support  vector  machine  to  fuel 
cell  control  system.  Schumacher  et  al.  [7]  proposed  a  water  man¬ 
agement  of  PEMFC  method  using  fuzzy  control.  According  to  the 
experimental  data,  Almeda  and  Smoesm  [8]  proposed  an  artificial 
neural  network  control  method  to  control  output  voltage  of  fuel  cell 
and  optimize  the  parameters  in  the  system.  Pukrushpan  co-workers 
[9,10]  used  feed-forward  and  feedback  strategies  to  control  the  flow 
rate  of  compressor  in  the  PEMFC  air  supply  system.  However,  the 
existing  control  approaches  used  for  PEMFC  were  based  on  linear 
models  which  were  linearized  at  a  specific  operating  point.  Chiu  et 
al.  [11]  proposed  a  linear  PEMFC  models  which  was  used  Jacobian 
linearization  via  a  Taylor  series  expansion  at  the  nominal  operating 
point.  However,  the  proposed  model  could  not  easily  achieve  sat¬ 
isfactory  dynamic  performance  under  large  disturbances  because 
of  the  operational  parametric  uncertainties  such  as  the  uncertain¬ 
ties  of  parametric  coefficients  for  each  cell  on  kinetic,  thermody¬ 
namic  and  electrochemical  foundations,  and  the  resistivity  of  the 
membrane  for  the  electron  flow.  Na  co-workers  [12,13]  presented 
a  nonlinear  controller  which  was  designed  based  on  the  nonlinear 
model  to  prolong  the  stack  life  of  fuel  cells.  The  simulation  results 
showed  that  the  proposed  nonlinear  controls  had  better  transient 
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performances  than  the  linear  controls.  However,  the  feedback  lin¬ 
earization  approach  was  used  to  only  achieve  the  local  linearization 
without  the  nonlinear  coordinate  transformation  because  of  lacking 
the  relative  degree  of  the  system.  The  proposed  nonlinear  con¬ 
troller  could  not  guarantee  the  robustness  with  the  operational 
parametric  uncertainties  and  the  internal  dynamics  problem  would 
be  appeared  out.  Therefore,  an  accurate  nonlinear  dynamic  model 
needs  to  be  developed  for  PEMFC  as  well  as  an  advanced  controller 
design  technique,  considering  the  nonlinearity  and  uncertainty  that 
need  to  be  proposed. 

Recently,  the  state  feedback  exact  linearization  for  nonlinear 
dynamic  models,  a  well-known  nonlinear  approach,  has  been 
widely  used  to  achieve  more  robust  transient  behavior  [14,16,17]. 
The  state  feedback  exact  linearization  uses  a  nonlinear  transforma¬ 
tion  to  transform  an  original  nonlinear  dynamic  model  into  a  linear 
model  by  diffeomorphism  mapping.  An  optimal  control  theory  is 
also  applied  to  obtain  a  linear  control  that  is  transformed  back  to  the 
original  space  by  using  the  nonlinear  mapping.  The  purpose  of  the 
dynamic  extension  algorithm  is  to  construct  an  extended  system  in 
which  the  relative  degree  is  equal  to  the  order  of  the  system  [15]. 
After  using  the  dynamic  extension  algorithm,  the  internal  dynamics 
problem  can  be  avoided  and  the  nonlinear  system  becomes  global 
linearized  via  the  state  feedback  exact  linearization. 

Hoc  control  approach  has  a  broad  base  of  support  because  of  its 
robustness  to  uncertainties  and  reliable  design  algorithms  [18-20]. 
The  weighting  functions  chosen  to  shape  the  sensitivity  func¬ 
tions  are  obtained  through  analysis  of  the  uncertainties  present 
in  the  system  as  well  as  from  frequency-domain  and  time-domain 
requirements.  Therefore,  in  the  framework  of  Ho 0  mixed  sensitivity 
design,  the  nonlinear  Hoo  robust  control  based  on  the  state  feedback 
exact  linearization  could  make  the  system  possess  better  robust 
stability. 

In  this  paper,  a  MIMO  dynamic  nonlinear  model  that  is  appro¬ 
priate  for  developing  a  nonlinear  robust  controller  is  proposed.  The 
state  feedback  exact  linearization  with  the  dynamic  extension  algo¬ 
rithm  is  applied  to  design  the  proposed  robust  controller,  based 
directly  on  the  nonlinear  dynamic  PEMFC  model.  The  dynamic 
response  of  the  proposed  model  is  tested  by  fuel  cell  test  system. 
The  control  law  obtained  from  the  state  feedback  exact  lineariza¬ 
tion  is  expected  to  be  more  robust  in  the  presence  of  large  dis¬ 
turbances  when  the  vehicle  running  process  is  simulated.  Further¬ 
more,  PEMFC  life  can  be  prolonged  and  protected  by  minimizing 
the  deviations  between  the  hydrogen  and  oxygen  partial  pressures. 


2.  The  modeling  of  PEMFC 


e~  Load  e~ 
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Fig.  1.  The  work  process  and  reaction  principle  of  the  PEM  fuel  cell. 


(3)  Due  to  a  slow  response  time  regarding  the  stack  temperature, 
the  operating  stack  temperature  is  assumed  to  be  constant. 

(4)  The  temperature  of  hydrogen  inside  the  anode  and  oxygen 
inside  the  cathode  are  equal  to  the  stack  temperature. 

( 5 )  For  water  management,  it  is  assumed  that  the  liquid  water  evap¬ 
orates  into  the  cathode  or  anode  gas  if  humidity  on  either  side 
drops  below  100%. 

(6)  It  is  assumed  that  pure  hydrogen  (99.99%)  is  fed  to  the  anode, 
and  air  that  is  uniformly  mixed  with  nitrogen  and  oxygen  by  a 
ratio  of  21:79  is  supplied  to  the  cathode. 

2.2.  The  output  voltage  model 

The  ideal  standard  potential  of  a  hydrogen/oxygen  fuel  cell 
which  is  obtained  for  standard  state  conditions  (25  °C  and  1  atm)  is 
1.229  V  with  liquid  water  product.  However,  the  actual  cell  potential 
is  decreased  from  its  equilibrium  potential  because  of  irreversible 
losses.  There  are  three  types  of  the  irreversible  losses,  namely  acti¬ 
vation  losses,  ohmic  losses  and  concentration  losses.  At  low  current 
densities,  the  activation  losses  of  oxygen  reduction  reaction  are 
almost  entirely  responsible  for  the  potential  drop  of  the  cell;  at 
the  high  current  densities,  the  concentration  losses  become  more 
significant.  The  output  voltage  of  the  single  cell  is  given  by  Eq.  (1) 
according  to  the  PEMFC  output  characteristics  empirical  equation 
which  is  developed  by  Amphlett  et  al.  [24]  and  Kim  et  al.  [25], 


2.1.  The  theory  of  PEMFC 


Uce\\  —  FNernst  —  ^act  —  ^ohmic  —  ^con 


(1) 


The  fuel  cell  is  an  electrochemical  energy  device  that  converts 
the  chemical  energy  of  fuel  directly  into  electricity  and  heat  with 
water  as  a  byproduct  of  the  reaction.  The  hydrogen  and  oxygen 
work  as  fuel  and  oxidant  and  they  need  to  be  humidified  before 
they  are  fed  into  the  cell/stack.  The  positively  charged  protons  dif¬ 
fuse  from  anode  through  one  side  of  the  membrane  and  migrate 
toward  the  cathode  [21-24].  The  electrons  pass  from  the  anode 
to  the  cathode  through  an  exterior  circuit  and  provide  electric 
power  along  the  way;  the  byproduct  water  will  be  drained  out  with 
the  gases.  The  work  process  and  reaction  principle  are  shown  in 
Fig.  1. 

In  order  to  simplify  the  analysis,  several  assumptions  are  made 
as  below: 


(1 )  The  stack  is  well  designed  so  that  all  cells  perform  similarly  and 
can  be  lumped  as  a  stack. 

(2)  All  gases  obey  the  ideal  gas  law  and  are  equably  distributed. 


where  FNernst  is  thermodynamic  potential,  Hact  is  activation  losses, 
Hohmic  is  ohmic  losses,  HCOn  is  concentration  loss. 


FNernst  — 


1 

2 F 


AG  -  A 5(1  -  Tref)  +  RT(lnPH2  + 


InPoj 
2  ‘ 


(2) 


where  AG  is  Gibb’s  free  energy  change,  AS  is  standard  mole  entropy 
change,  R  is  gaseous  constant,  F  is  Faraday’s  constant,  Tref  is  refer¬ 
ence  temperature,  Ph2  is  partial  pressure  of  hydrogen,  Po2  is  partial 
pressure  of  oxygen. 

The  activation  loss  of  PEMFC  is  caused  by  the  sluggish  kinetics  of 
the  reactions  taking  place  on  the  active  surface  of  electrodes  [26,28]. 

Ua  ct  =  +  H2T  +  ^F(lnCO)  +  §4T(ln/s)  (3) 

where  Jst  is  the  load  current  of  PEMFC,  £i,  §2.  £3.  §4  are  the  model 
coefficient  got  by  experimental  data  fitting  based  on  electrochem¬ 
istry,  thermodynamics  and  fluid  mechanics,  CO2  is  dissolved  oxygen 
concentration  in  the  interface  of  the  cathode  catalyst. 
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According  to  the  definition  of  Henry’s  law  [27],  it  yields, 


5.08  x  106exp(— 498/T)  '  ' 

Different  from  the  ohmic  losses  empirical  equation  which  is 
already  presented,  in  this  model  ohmic  losses  consist  of  the  voltage 
drop  which  is  caused  by  Pm,  the  equivalent  membrane  impedance 
and  the  voltage  drop  that  is  caused  by  Rc,  the  contact  resistances 
both  between  the  membrane  and  electrodes  as  well  as  the  elec¬ 
trodes  and  the  bipolar  plates  which  is  constant  once  the  cell  is 
fabricated  [27,28]. 


where  Ph2  and  Pw, an  are  the  partial  pressure  of  hydrogen  and  water 
in  the  anode  respectively,  WHl  in  is  the  inlet  mass  flow  rate  of  hydro¬ 
gen,  Wh2 .purge  is  the  hydrogen  mass  flow  rate  purged  out  the  anode, 
WH2,ret  is  the  reacted  mass  flow  rate  of  hydrogen,  Wvanin  is  the 
water  mass  flow  rate  into  the  anode,  Wv,an, purge  is  the  purged  water 
mass  flow  rate,  Wvmbr  is  the  water  flow  rate  across  the  membrane. 

With  the  measured  inlet  flow  rates  and  the  stack  current,  the 
outlet  flow  rates  are  given  by  Ain  which  is  the  summations  of  anode 
inlet  flow  rates  minus  the  reacted  gas. 

WVet  =  N^=Qjist  (14) 


^ohmic  —  htiR-M  +  ^c) 


(5)  Wh2  .purge  —  (An  —  Q-\ht)pH2 


(15) 


(6) 


where  l  is  the  thickness  of  the  membrane,  Afc  is  the  cell  active  area. 
The  resistivity  rM  of  Nafion  series  proton  exchanging  membrane  can 
be  calculated  by  Eq.  (7)  [27]. 


^H2,in  —  ^H2An 
Ph2 


%  = 


Ph2  +  Pw,a 
Ain=SH2Uj<a 


(16) 

(17) 

(18) 


=  181.6[l+0.03(/st/afc)  +  0.062(T/303)2(/stMfc)2'5] 

M  [A.  -  0.634  -  3(/st/.Afc)]exp[4.18(T  _  303/7)]  1  ’ 

where  X  is  the  water  content  of  the  membrane. 

Concentration  losses  are  caused  by  the  mass  transportation 
which  in  turn  affects  the  concentration  of  the  hydrogen  and  oxygen 
at  high  current  density  [28,30]. 

Ucon=Bln(\-fL.)  (8) 

where  B  is  determined  by  the  PEMFC  and  its  working  status,  ist  is 
the  current  density  of  fuel  cell,  imax  is  maximum  current  density. 

There  exists  charge  double  layer  in  the  PEMFC,  Larminie  et  al. 
[1]  and  Wang  et  al.  [30]  analyzed  the  charge  double  layer  at  the 
surface  of  the  cathode.  The  equivalent  capacitance  C  can  smooth 
the  voltage  drop  effectively.  As  the  function  of  the  charge  double 
layer,  the  PEMFC  bears  eminent  dynamic  characteristics.  If  ud  is  the 
overall  voltage  drop,  the  dynamic  characteristics  of  the  single  cell 
can  be  expressed  in  the  differential  Eq.  (9),  which  is  distinct  from 
the  steady-state  model. 


where  N  is  the  cell  number  Yh2  is  the  mole  fractions  of  hydrogen, 
Fh2  is  the  pressure  fraction  of  hydrogen  inside  the  anode,  ua  is  the 
input  control  variable  of  anode,  ka  is  the  conversion  factors  of  anode, 
Sh2  is  hydrogen  stoichiometry. 

The  membrane  hydration  model  captures  the  effect  of  water 
transport  across  the  membrane.  Both  water  content  and  mass  flow 
are  assumed  to  be  uniform  over  the  surface  area  of  the  membrane 
[31,32].  The  water  inlet  flow  rates  on  the  anode  and  the  cathode  are 
expressed  in  terms  of  the  relative  humidity,  saturation  pressure, 
and  total  pressure  on  each  side. 


i^v.mbr  —  NMvAfcUd  p  —  Q^st 

(19) 

Wv,an, purge  —  (An  —  Q2ist)Fv,an 

(20) 

W  —  a 

v’an'ln  PH2+Pw,an-<PaPv  ln 

(21) 

P  Pw,an 

V’an  “Pu  +P 

1  H 2  '  1  w,an 

(22) 

duo  _  ht  _  Ud  (Qs 

dt  ~  C  x  1  j 

where  the  time  constant  r  will  change  with  the  load  changing, 
which  will  affect  the  dynamic  response  of  the  voltage  as  shown 
in  Eq.  (10). 

x  =  e(Uact/+tUcon)  (10) 

Considering  synthetically  the  thermal  characteristic,  momen¬ 
tum  characteristic  and  the  mass  transportation,  the  output  voltage 
of  PEMFC  can  be  expressed  as  Eq.  (11). 

t^cell  =  pNersnt  —  t/ohmic  —  (H) 


2.3.  The  anode  pressure  model 


According  to  the  ideal  gas  law  and  the  mode  conservation  rule, 
the  partial  pressure  of  each  gas  is  balanced  by  the  gas  inlet  flow  rate 
minus  the  gas  consumption  and  the  gas  outlet  flow  rate.  The  partial 
pressure  derivatives  are  given  as  follows  [29,31,32]. 

Anode  mole  conservation: 


dP»2 

dt 

dPw.an 

dt 


Tl_(^H2,in  —  Wh2  .purge  —  k^H2,ret) 

va 

RT 

-  y-(Wv,an,in  -  Wv,an, purge  -  Wv.mbr) 


(12) 

(13) 


logio(Pv)  =  -1.69  x  10-10T4  +  3.85  x  lCT7!3  -3.39  x  lCT4!2 

+  0.143T- 20.92  (23) 

where  nd  is  the  electro-osmotic  drag  coefficient,  Mv  is  the  vapor 
molar  mass,  (pa  is  the  relative  humidity  in  the  anode,  Pv  is  the 
saturation  pressure. 


2.4.  The  cathode  pressure  model 


This  model  captures  the  cathode  air  flow  behavior.  Cathode  mole 
conservation  [29,31,32]: 


=  ^(W°2.in  “  W02,out  -  W°2,ret) 

^  =  £(wN2,in  -  wN2,0Ut) 

— ^ —  =  ^“(Wv.ca.in  —  Wv,Ca,out  +  Wv.gen  +  Wvmbr) 


(24) 

(25) 

(26) 


where  P02 ,  Pn2  and  PW)C a  are  the  partial  pressure  of  oxygen,  nitrogen 
and  water  in  the  cathode  respectively,  W02.in  and  WN2  in  are  the 
oxygen  mass  flow  rate  and  the  nitrogen  mass  flow  rate  into  the 
cathode  respectively,  W02,ret  is  the  reacted  oxygen  mass  flow  rate  in 
the  cathode,  Wo2,0ut  and  WN2i0Ut  are  the  oxygen  mass  flow  rate  and 
the  nitrogen  mass  flow  rate  out  the  cathode  respectively,  Wv,gen  is 
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the  reacted  water  mass  flow  rate,  WViCa>in  and  WV)Ca,out  are  the  water 
mass  flow  rate  into  the  cathode  and  out  the  cathode  , respectively. 


W  MAfc’SC  i  Ql 

Wo2,ret  =N  4p  =  >st  2 

(27) 

Wo2, out  =  ^Qn  -  F02 

(28) 

l^N2,out  —  Qn^N2 

(29) 

^02,in  —  ^02  Qn 

(30) 

^Nj.in  =  ^-in 

(31) 

F  P(h 

2  Po2  +  Pn2  +  Pw,ca 

(32) 

Pn9 

n  _  1N2 

2  Po2  +  Pn2  ~\~Pw,ca 

(33) 

Qn  —  5airUc/<c 

(34) 

where  Cin  is  the  summations  of  cathode  inlet  flow  rates,  Yq2  and 
Yn2  are  the  mole  fractions  of  oxygen  and  nitrogen,  Fq2  and  FN2  are 
the  pressure  fraction  of  oxygen  and  nitrogen  inside  the  cathode,  uc 
is  the  input  control  variable  of  cathode,  kc  is  the  conversion  factors 
of  cathode,  sair  is  air  stoichiometry. 

Wv.ca.out  —  (Qn  —  Q2ht)Fv,ca 

(35) 

\a/  <PcPv  r 

v,ca,‘n  _  po,  _|_  pw  ca  —  y>cPv  m 

(36) 

defined  input  v.  The  notation  and  concepts  of  differential  geom¬ 
etry  are  essential  to  understand  this  approach.  Consider  a  MIMO 
nonlinear  system 

m 

*=/(*> +5><*>*«,  (39) 

i=  1 

yi  =  hi{x),  i  =  l,2,  ...,m 


where  x  is  an  n-dimensional  state  vector  that  is  assumed  to  be 
measurable,  u  and  y  are  m-dimensional  input  and  output  vectors. 

The  Lie  derivative  of  a  scalar  function  h(x)  with  respect  to  a  vector 
function  f[x)  is  defined  as 

Lfh(x)  =  Vhf  =  ^/(x)  (40) 

Repeated  Lie  derivatives  can  be  defined  recursively  as 
LfK*)  =  h)  =  V(L)-1  h)f 

L°h(x)  =  h(x)  k  =  1,2,...  1  J 

Similarly,  in  the  case  of  another  vector  field  g 


LgLfh(x)  =  V(Lfh)g  (42) 

The  output  needs  to  be  differentiated  for  r  times  until  it  is 
directly  related  to  the  input  u.  The  number  r  is  called  the  relative 
degree  of  the  system.  The  MIMO  system  is  said  to  have  a  vector 
relative  degree  {r\ ,  r2, . . .  ,rm}  at  a  point  x0  if 


Wv,gen 
Tv.ca  = 


=  N^=(hUt 

Pw,ca 

Pq2  +  f*N2  +  ^w,ca 


(37)  (1 )  LgjL^hiix)  =  0  for  all  1  <j  <  m ,  all  1  <  i<  m,  k<  rx  - 1,  and  for  all  x 

in  the  neighborhood  of  x0; 

(38)  (2)  m  x  m  matrix 


where  (pc  is  the  relative  humidity  in  the  cathode. 

In  addition,  if  the  mass  of  water  calculated  is  greater  than  the 
maximum  mass  of  vapor,  the  liquid  water  formation  occurs  simul¬ 
taneously.  The  flow  rates  of  liquid  water  leaving  the  anode  and 
cathode  are  dependent  upon  the  saturation  state  of  each  gas  [10]. 
The  mass  of  liquid  water  and  vapor  is  calculated  as  follows: 

ifmWianorca  <  tnv,max,an  or  ca  tnv,anorca  —  tnw.anorca?  ^l.anorca  — 

ifmWianorca  >  nTv,max  an  or  ca  tnv,anorca  —  tttw,max,an  or  ca  > 

^l.anorca  —  tMw>anorca  —  tnw,max,anorca- 

In  designing  PEMFC  control,  the  main  focus  is  on  the  control  of 
hydrogen  and  oxygen  partial  pressures,  which  can  avoid  unwanted 
pressure  fluctuation  and  prevent  the  membrane  electrode  assem¬ 
blies  (MEAs)  from  collapsing  by  minimizing  the  pressure  difference 
between  the  anode  and  the  cathode. 

3.  Nonlinear  control  by  state  feedback  exact  linearization 

For  chemical  process  control,  nonlinear  control  theory  devel¬ 
oped  from  differential  geometry,  known  as  exact  linearization  or 
feedback  linearization,  has  more  attractive  because  many  chemical 
processes  are  basically  of  high  nonlinearity  [31,33].  Hence,  one  of 
the  main  motivations  of  utilizing  state  feedback  exact  linearization 
for  a  PEMFC  system  is  inherently  a  nonlinear  chemical  process. 

3.L  State  feedback  exact  linearization 


■  LgjLp-’hjCx) 


E(x)  = 


Lg^-’Mx) 


LgmLp-'h lM  ' 


LgmLjm  1  hm(x) 


(43) 


is  nonsingular  at  x  =  x0,  which  is  called  as  a  decoupling  matrix.  If 
E(x)  is  singular,  a  considerably  more  elaborate  analysis  is  required. 

Considering  a  nonlinear  system  has  the  relative  degree 
V\  +  *2  +  •  •  •  +  rm  =  n,  a  nonlinear  coordinate  transformation  in  the 
state  space  is  defined  as 


z[  =  0\x)  = 

1 

2  2  • 

■ 

i _ 

_ 

h,(x) 

Lfhj(x) 

yk(x)_ 

_Lkf-'hi(x)_ 

1  <  i  <  m,  1  <  k  <  Tj  (44) 


where  is  called  a  local  diffeomorphism,  in  which  the  map  between 
the  new  input  v  and  the  output  is  exactly  linear  for  all  x  in  the 
neighborhood  of  x0.  Then  there  exists  a  nonlinear  static  feedback 
such  that  the  closed-loop  system  in  the  new  coordinates  is  linear 
and  controllable. 

In  the  new  coordinates  the  system  is  given 
z\  =  z\ ,  1  <  i  <  m 

by  =zv  1  <i<m  where 

m 

z'r,  =  bj(z)  +  y>(z)u,,  1  <  i  <  m 

j=i 


The  objective  of  state  feedback  exact  linearization  is  to  create 
a  linear  differential  relation  between  the  output  y  and  a  newly 


b,(z)  =  Lp/1,(0 -1  (z))  ay(z)  =  LgjLp-i  hi :(</>-’  (z)) 
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Suppose  that  the  matrix  z\  (1  <  i  <  m)  is 


’Zr,  " 

'bj(z)' 

=  b(z)  +  E(z)u  = 

3Z. 

_bm(z)_ 

'an(z)  ■■■  aim(z)" 

_am  l(z)  Qmm(z)_ 

_  _ 

Then  a  newly  vector  input  is  defined  based  on  above  the  matrix 


that  the  nonlinear  system  is  global  linearized  and  has  a  satisfactory 
controller. 

The  purpose  of  the  dynamic  extension  algorithm  is  to  construct, 
starting  with  a  system  in  which  the  relative  degree  r  =  r\  +  r2. . .  +  rm 
is  not  equal  to  n,  an  extended  (and  feedback-modified)  system  in 
which  the  relative  degree  is  possibly  larger,  and  therefore-possibly 
after  a  number  of  iterations  -  equal  to  n. 

It  is  possible  to  find  two  integers  i0,  j0,  and  make  an  element  of 
E(x)  that  is 

aia)0(x°)  =  Lgj0Lpo-'  hjg(x0)  ±  °  (51) 

To  define  the  dynamic  feedback,  in  which  p(x)  and  q(x)  are  arbi¬ 
trary  function  satisfying  p(x°) =  0  and  q(x°)  =  1. 


b(z)  +  £(z)u  =  v  = 


"y(!ri)" 

Vm_ 

- 1 

••  3e 

_ i 

(47) 


Therefore  the  state  feedback  control  law  of  the  MIMO  nonlinear 
system  is  defined  as 

u  =  -£_1  (z)b(z)  +  £_1  (z)  v  (48) 


Note  that  the  control  law  in  Eq.  (48)  transforms  the  nonlinear 
system  into  a  linear  one  in  which  the  aforementioned  input-output 
relation  is  linearized  and  decoupled. 

The  resulting  closed-loop  system  which  is  so-called  Brunovsky 
canonical  form  is  then  governed  by  the  equations 

z\  =  z\,  1  <  i  <  m 

: .  (49) 

z\  1  =  zl. ,  1  <  i  <  m 

ri~ 1  ri  ~  ~ 

z\.  =  Vj,  1  <  i  <  m 


Uj=Vj,  j^j0 

1  m 

Ui°  =  %~M{p{-x)+q{x^~  53  a<dW 


j= i 
j  *}  0 


(52) 


The  composition  of  Eqs.  (39)  and  (52)  define  a  new  system  that 


(53) 


III  f  V  III 

X  =/(*)+  V'  gj(.x)Vj  +  gj  7. (p(x)  +  q(x)g -  V  aiol(x)Vj] 
aia/oW 

j  =  1  J  =  1 

j  =/=jo  j  ¥=j0 

s  =  vk 

yt  =  h(x),  i  =  1,  2, . . . ,  m 


By  definition,  the  rioth  derivative  of  yi0  can  be  expressed  in  the 
form 


Above  the  coordinates  transformation  and  state  feedback  prob- 
lems  are  so-called  state  feedback  exact  linearization  problem.  Then  3%'°  =  ^/l°  hi0  M  +  /  Vq/Qj(x)uj 
the  new  linear  state-space  equation  is  j= 1 


(54) 


z\.  =  AjZ}..  +  Bji/j,  1  <  i  <  m  (50) 

with 


"0 

1 

0  ••• 

0" 

"0" 

0 

0 

1  ••• 

0 

0 

Bi  = 

0 

0 

0 

0  ••• 

1 

_0 

0 

0  0 

0_ 

mxm 

_1_ 

The  matrix  Az  and  the  vector  Bz  satisfy  the  controllable  condition. 
rank(Bz  AtB \ . . .  A"-1^,-)  =  n 

In  this  paper,  the  state  feedback  exact  linearization  is  utilized 
according  to  the  MIMO  dynamic  nonlinear  model  of  the  PEMFC. 

3.2.  Dynamic  extension  algorithm 

To  apply  the  state  feedback  exact  linearization,  the  relative 
degree  vector  has  to  be  checked  such  that  the  internal  dynamics 
problem  is  avoided.  Without  solving  the  internal  dynamics  issue, 
the  nonlinear  controller  would  be  practically  meaningless  because 
the  instability  of  the  internal  dynamics  causes  undesirable  phe¬ 
nomena  such  as  a  poor  response  and  critical  damage  to  the  system. 
Therefore,  the  relative  degree  associated  with  the  state  feedback 
exact  linearization  must  be  the  same  as  the  order  of  the  system  so 


Thus,  choosing  p(x)  =  -Lp0  hio  (x)  and  q(x)  =  1  in  the  law  yields, 
for  the  rlQth  derivative  ofyl0,  the  simple  expression 

4r'o}=?  (55) 

The  latter  in  turn  yields 

y^o+D  =  vj0  (56) 

Therefore,  in  the  composed  system  Eq.  (53),  the  lowest  deriva¬ 
tive  of  yiQ  which  explicitly  depends  on  the  input  is  precisely  the 
(rlQ  +  l)th  derivative. 

Suppose  the  dynamic  extension  algorithm  has  been  iterated  k 
times.  Let 

X  =f{x)  +  g(x)[H(x,  |)  +  K(x,  £)P] 

£  =  F(x,  I)  +  G(x,  £)t> 

(57) 

v  =  a(x,  $)  +  fi(x,  f)v 
y  =  h(x) 

with  f  e  Rk,  denote  the  composition  of  the  fcth  feedback  laws  of  the 
form  constructed  at  each  stage  of  the  algorithm. 

After  iterating  k  times  of  the  dynamic  extension  algorithm,  the 
relative  degree  r  is  equal  to  n.  Therefore,  there  is  no  internal  dynamic 
problem,  and  the  nonlinear  system  becomes  global  linearized  via 
the  state  feedback  exact  linearization. 
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4.  Nonlinear  robust  control  of  PEMFC 

In  this  section,  a  MIMO  dynamic  nonlinear  model  of  PEMFC  is 
developed  by  Eqs.  (l)-(38),  and  it  is  then  used  to  design  a  nonlin¬ 
ear  robust  controller  by  adopting  state  feedback  exact  linearization 
in  order  to  minimize  the  pressure  deviation  between  the  hydrogen 
and  oxygen.  The  main  purpose  of  keeping  the  difference  of  pres¬ 
sures  in  a  certain  small  range  is  to  avoid  a  detrimental  degradation 
of  the  fuel  cell  voltage,  because  the  fuel  cell  voltage  is  a  function  of 
the  pressures,  each  pressure  needs  to  be  appropriately  controlled, 
and  the  pressures  have  bigger  impact  on  the  performance  of  fuel 
cells  than  other  parameters  [9,13].  In  addition,  keeping  the  small 
difference  of  pressures  is  also  to  protect  the  membrane  from  dam¬ 
age,  and  therefore,  prolong  the  fuel  cell  stack  life.  In  this  paper,  a 
pressure-control  optimization  algorithm  is  focused  on  developing 
for  the  PEMFC  system.  The  stack  current  is  considered  as  a  distur¬ 
bance  to  the  system  instead  of  an  external  input. 

4.1  The  state-space  MIMO  dynamic  model 

Consider  the  following  MIMO  nonlinear  system  with  a  distur¬ 
bance: 


m 

x=f(x)  +  'Y)gi(x)ui  +  L(x)w 
1=1 

yi  =  hi{x),  i  =  1, . . .  2,  m 

where  xcRn  is  the  state  vector,  UcRm  is  the  input  vector,  ycRm  is 
the  output  vector,  and  /(x)  and  g(x)  are  n-dimensional  smooth  vec¬ 
tor  fields.  The  w  represents  the  disturbance  variables,  and  L(x)  the 
dimensional  vector  field  directly  related  to  the  disturbance. 

According  to  Eqs.  (l)-(38),  the  state-space  MIMO  nonlinear 
dynamic  system  model  of  PEMFC  is  expressed  as  follows: 


L(x)  =  RT 


Qi  Qi*i 
Va  +  Va(*l  +  *2) 
02*2  02 
Va(*l+*2)“  Va 
Q.1  0.1*3 

2Vc  2Vc(x3+X4  +  x5) 
0 


Ol _ 0.1*5 _ 02*5 _ ^  Q2 

-  Vc  ^c(*3  +  *4  +  *5)  VC(X3  +  *4  +  *5  J  Vc  _ 


Eq.  (59)  implies  that  the  input-output  behavior  of  the  PEMFC 
system  is  nonlinear  and  coupled.  Each  control  variable  u  shows 
up  after  the  first  derivative  of  each  y ^  =X\  and  y2  =*3,  the  relative 
degree  vector  {r\ ,  r2}  is  {1, 1 },  so  r=  r\  +  r2  =2  <  5,  and  the  decoupling 
matrix  E(x)  is  defined  as 


'Lglhi(x)  Ig2hi(x)' 
_Lglh2(x)  Lgl  h2{x)  _ 


(60) 


In  order  to  achieve  the  control  objective,  three  steps  need  to 
follow. 


( 1 )  Applying  the  dynamic  extension  algorithm,  an  extended  system 
in  which  the  relative  degree  r  is  equal  to  5  after  a  number  of 
iterations  is  achieved  such  that  the  internal  dynamics  problem 
is  avoided. 

(2)  Using  a  nonlinear  coordinate  transformation  in  the  state-space, 
there  exists  by  a  nonlinear  static  feedback  such  that  the  closed- 
loop  system  can  be  changed  into  Brunovsky  canonical  form  in 
which  the  system  is  linear  and  controllable. 

(3)  Obtaining  a  nonlinear  robust  control  law  that  not  only  can  com¬ 
pensate  nonlinearities  but  also  can  decouple  and  linearize  the 
input  and  output  behaviors. 


X=/(x)+g!(x)U!  +  g2{x)u2  +  L(x)W 


>r 

"*r 

'hiW 

J2. 

.*3. 

MW. 

where 


*r 

-P»2  - 

*2 

fw,an 

*3 

= 

Po2 

*4 

Pn2 

*5- 

-  Pw,c  a  - 

y\ 

Ph2 

y2  _ 

lP°2\ 

Ua 

u2 

Uc 

w  =  im,  /(x)  =  0 


gi(x)  =  K7sH2 


Yn2ka  _  /<a  *1 

Va  Va  Xi  +  X2 
(PaPvka  /<a  *1 

Va(Xi  +  *2  —  ^a^v)  Va  *1  +  *2  ’ 

0 

0 

0 


&(x)  =  RTsa[] 


0 

0 


Yo2  1<c  kc  x3 

Vc  Vc  x3  +  x4  +  X5 

Yn2  kc  kc  x4 

Vc  Vc  X3  T  X4  T  X5 

(PcPykc  kc  X5 

Vc(x3  +  X4  +  X5  -  (pcPv )  Vc  X3  +  X4  +  X5 


(59) 


/  lCCOlUlll^  LU  I— f  U1I 

extension  algorithm,  the  extended  system 


is  as  roiiows, 


x  = 


=  /(x)+gi(x)v1  +  g2(x)v2  +  L(x)w 


"yr 

■*r 

'hi(x)' 

_y2. 

_*3  _ 

MW. 

where 


/(*)= 


S12 1 
§3 

624  t 

^ 3 

625  t 

§2 

0 

0 


>gl(*)  = 


"0" 

"0" 

0 

0 

,Mx)  = 

1 

0 

_0_ 

_1_ 

(61) 


(62) 


After  the  iteration,  the  relative  degree  vector  f  is  equal  to  5. 
Then  the  decoupling  matrix  F(x)  and  the  nonlinear  control  law  are 
defined  as 


Lg1Lfh-i(x)  Lg2Ljrh-i(x) 

Lgi  Lfh2(x)  TglT/h2(x) 


E\x)  = 


(62) 
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Fig.  2.  Structure  of  mixed  sensitivity  problem. 


Ul 

U2 


(l/a11(x))(p1(x)  +  q,(x) 


(l/a22(x))(p2(x)  +  (}2(x)  /  v2) 


I  (l/a'^x))  J 

/ ' 


V\) 


-E~\x) 


Ln(x) 

MsW 


w 


(63) 


Therefore,  there  is  no  internal  dynamic  problem,  and  the  nonlin¬ 
ear  system  becomes  global  linearized  via  the  state  feedback  exact 
linearization.  The  nonlinear  coordinate  transformation  in  the  state- 
space 


Xl" 

ZJ  = 

= 

_03M_ 

_§2_ 

"*3" 

_?3  . 

(64) 


Then  the  Brunovsky  canonical  form  is  obtained 


"4 

’4’ 

z1  =  Al 

Z2 

A. 

+  B1V1  yi=C! 

4 

.4. 

z2  = 
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"4" 

A\ 
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1 

0" 
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"r 
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0 

0 

1 

,  Bi  = 
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0 

4.2.  Hoo  robust  optimal  control 


(65) 


To  begin  with,  consider  a  standard  mixed  sensitivity  problem 
[20]  shown  in  Fig.  2.  Where  r,  e,  u,  d  and  y  are  the  reference  input 
signal,  the  tracking  error  signal,  the  control  input  signal,  the  dis¬ 
turbance  signal  and  the  feedback  output  signal,  respectively.  I<{s)  is 
the  controller,  and  G(s)  the  plant. ,y2, y3  is  the  weighting  estimate 
signal  of  generalized  plant. 

Using  the  frequency-dependent  weighting  functions  Wi(s), 
W2  (s)  and  W3  (s)  to  weight  the  signals  e,  u  and y  respectively,  a  gener¬ 
alized  system  is  constructed,  which  can  be  furthermore  converted 
into  the  standard  small  gain  problem.  The  mixed  sensitivity 


design  problem  is  to  find  a  proper  rational  controller  that  stabilizes 
the  closed-loop  system  and  satisfies  [20]: 


min  0 


WiS 

W2R 

w3t 


=  Y 


00 


(66) 


where  y  is  the  optimal  H ^  norm,  W\  (s),  W2(s)  and  W3(s)  are  weight¬ 
ing  functions  for  the  sensitivity  matrix  function  S(s),  the  control 
signal  sensitivity  matrix  R( s),  and  the  complementary  sensitivity 
matrix  function  T(s),  respectively,  are  expressed  as: 


S  =  (/  +  G/C)-1 

R  =  K(I  +  GK)-1  =  KS  (67) 

T  =  GK(I  +  GK)"1  =  /  -  S 


For  Wi(s)  and  W3(s),  they  need  to  satisfy  the  inequality  con¬ 
straint: 

o'[W1“10'w)] +  o-[W3“1(jw)]  >  1  (68) 


As  it  is  known,  the  sensitivity  function  is  especially  interesting 
from  the  control  point  of  view  as  it  has  properties  that  charac¬ 
terize  both  the  quality  and  the  speed  of  the  closed-loop  temporal 
response.  Bearing  in  mind  that  Wi(s)  must  weight  the  sensitivity 
function,  its  design  is  proposed  to  be  a  diagonal  matrix  of  transfer 
functions.  And  each  diagonal  element  of  Wi(s)  must  be  designed 
bearing  in  mind  that  its  inverse  should  shape  to  the  each  diagonal 
element  of  5(s)  as  an  upper  bound.  The  inverse  of  Wi  (s)  must  be  sta¬ 
ble,  minimum  phase  and  the  maximum  singular  value  of  yW\ -1  (s) 
must  greater  than  the  maximum  singular  value  of  S(s)  for  overall 
frequency  range: 

cr[S(jw)]  <  cr[yW^{jw)]  (69) 

Since  the  class  of  uncertainty  employed  corresponds  to  the  mul¬ 
tiplicative  output  one,  W3(s)  associate  robust  stability  condition 
and  the  maximum  singular  value  of  yW3_1(s)  must  greater  than 
the  maximum  singular  value  of  T(s)  for  overall  frequency  range: 

a[T(jw)]  <a[yW3~'(jw)]  (70) 

Moreover,  taking  into  account  that  W3(s)  must  weight  to  the 
complementary  sensitivity  function,  it  is  desirable  that  the  module 
of  W3(s)  has  a  high  value  in  order  to  impose  that  T(s)  has  a  small 
gain  at  high  frequencies. 

W2(s)  is  usually  taken  as  diagonal  constant  matrix  to  avoid 
increasing  the  order  of  the  controller. 

According  to  the  requirements  on  the  robustness  and  perfor¬ 
mance  of  the  PEMFC  system  which  is  mentioned  above,  the  struc¬ 
tures  of  Wn,  Wi2,  Wi3,  W2i,  W22,  W23  are  selected: 


Wn(s) 

Wu(s) 

W13(s) 

W21(s) 

W22(s) 

W23(s) 


42.3  x(s  +  106.8) 
(s  +  0.0012)2 

0.0000188 
0.005(s  +  13.9) 

118.3  x(s  + 60.8) 
(S  +  1.15)2 

0.00021 
0.002(s+  10.6) 


(71) 


Meanwhile,  the  H ^  mixed  sensitivity  problem  could  be  tackled 
based  on  solving  2-Riccati  Eq.  (19).  Then  the  optimal  controllers 
for  inlet  flow  rates  of  anode  and  cathode  which  norm  y\  and  y2  are 
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Fig.  3.  Overall  control  scheme  of  PEMFC  system. 


0.85  and  0.58  are  obtained: 

v  ,  ,  2.8  x  106(s  +  0.46)(s  +  43.7)(s  +  2.2)(s  +  1.9) 

i\dS  J  =  - - - 

(s  +  3.2  X  104)(s  +  0.001  )(s  +  0.06)(s  +  101 ,3)(s  +  60.4) 

„  ,  .  11.2  x  105(s +  46.9)(s  + 0.01  )(s  + 0.2) 

I<2  (S)  =  - - - - - - — 

(S  +  2.384)(s  +  1 .9  x  104)(s  +  2.4)(s  +  1 .9) 

(72) 

Therefore,  according  to  Eq.  (61),  the  outputs  Ph2  and  Po2  track 
asymptotically  the  desired  trajectory  3  atm  by  adding  the  optimal 
Hoc  controllers  mentioned  in  [34,35].  As  shown  in  Fig.  3,  the  main 
objective  of  this  control  scheme  is  to  design  a  nonlinear  robust  con¬ 
troller  by  appropriately  defining  a  transformation  mapping  scheme 
that  transforms  the  original  nonlinear  system  into  a  linear  and  con¬ 
trollable  closed-loop  system. 

5.  Experimental 

The  PEMFC  test  system  was  set  up  in  the  Fuel  Cell  Application 
Centre  (FAC),  Temasek  Polytechnic  Engineering  School. 

5.1  The  Cell/Stack 

In  this  experiment,  a  20-cell  stack,  which  was  designed  and  fab¬ 
ricated  by  Singapore  GasHub  and  the  FAC,  as  shown  in  Fig.  4,  was 
tested  for  the  validation  of  simulation  results.  A  Nafion  112  mem¬ 
brane  was  used  as  the  electrolyte,  the  Pt  loading  on  electrodes  is 
0.4  mg  cm-2  in  the  catalyst  layer,  graphite  foil  with  serpentine  flow 
channel  was  used  as  the  bipolar  plate,  and  the  active  area  of  the 
electrode  is  150  cm2  for  the  stack  [36,37]. 

5.2.  Steady -state  polarization  measurement 

For  the  steady-state  polarization  measurement,  both  the  hydro¬ 
gen  and  air  are  humidified  before  entering  the  fuel  cell  stack.  The 
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Fig.  5.  Polarization  and  power  curve  of  the  experimental  PEMFC  stack. 


air  is  humidified  to  100%  by  using  EFH-100WA  solid-state  humid¬ 
ifier  from  EnerFuel,  Inc.,  while  a  hydrogen  recycling  pump  is  used 
to  achieve  the  self-humidification  of  hydrogen  at  the  anode,  as  pre¬ 
sented  in  Fig.  5.  The  stack  temperature  is  323  K,  the  pressure  of  both 
sides  is  0.5  atm,  and  the  fuel  cell  stack  repeatedly  runs  under  full 
load  to  achieve  stable  performance  before  the  steady-state  polar¬ 
ization  measurements  are  carried  out. 

5.3.  Transient  response  measurement 

To  study  the  transient  response  of  the  fuel  cell,  two  items  would 
be  clarified:  cold  start  and  hot  start.  In  this  experiment,  cold  start 
means  that  the  hydrogen  goes  to  the  ON/OFF  valve  only  at  the 
moment  the  system  starts  and  hot  start  means  that  the  hydro¬ 
gen  always  fully  fills  up  the  cell/stack  even  when  the  system  does 
not  work.  The  advantage  of  cold  start  is  that  the  gas  is  confined 
within  the  pipe  when  the  system  is  in  standby,  and  is  more  reli¬ 
able  and  safe;  the  disadvantage  is  that  it  will  affect  the  lag  time  of 
the  system.  The  advantage  of  hot  start  is  that  the  starting  of  the 
system  is  swift  and  the  lag  of  the  system  can  be  shortened  effec¬ 
tively;  the  disadvantage  is  that  the  fuel  stays  inside  the  cell/stack, 
and  some  reaction  and  crossover  between  the  MEA  will  take  place; 
meanwhile,  there  will  be  leakage  of  the  fuel,  and  thus,  the  system 
fabrication  must  be  designed  in  a  special  way.  If  a  small-capacity 
battery  or  an  ultra-capacitor  is  employed  in  the  output  port  of  the 
fuel  cell,  the  disadvantage  of  the  cold  start  can  be  compensated,  and 
it  is  accepted  by  the  user  mostly  for  its  safety.  Thus,  we  focus  on  the 
cold  start  in  this  study. 

5.4.  The  Equipments 

The  voltage  response,  the  voltage  versus  current  curve,  and  the 
resistance  of  the  cell/stack  are  measured  by  Scribner  Associates  890 
CL  and  LeCroy  Waverunner  LT344.  The  mass  flow  controllers  with 
electromagnetic  valves  are  installed  to  control  the  hydrogen  flow 
and  air  flow  to  record  the  response  of  the  output  voltage  and  mini¬ 
mize  the  delay  of  the  gas  pipe.  The  humidifiers  are  used  to  prevent 
dehydration  of  the  fuel  cell  membrane  and  improve  the  reaction 
rate.  A  hydrogen  purge  system  including  segregators  is  applied  to 
increase  the  utility  of  fuel.  The  water  cooling  system  used  is  Julabo 
FP40,  which  uses  PID  algorithm  to  control  the  temperature  and  has 
a  resolution  of  0.1  °C.  Fig.  6  shows  the  diagram  of  the  PEMFC  test 
system. 

6.  Model  validation  and  simulation  results 


Fig.  4.  Experimental  PEMFC  Stack. 


In  this  paper,  the  Matlab/Simulink  is  used  to  setup  the  PEMFC 
system  dynamic  model  with  nonlinear  robust  controller.  To  verify 
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Fig.  6.  Diagram  of  the  PEMFC  test  system. 


the  validity  of  model,  an  experimental  of  a  1  kW  PEMFC  stack  with 
20  cells  was  conducted.  The  performances  of  the  stack,  such  as  out¬ 
put  voltage  versus  current  ( V-I )  curve  and  power  versus  current 
(P-J)  curve,  are  presented  in  Fig.  5.  It  can  be  seen  that  the  stack  can 
give  a  1  kW  output  at  about  70  A  of  current. 

Fig.  7  shows  the  measurement  results  of  the  transient  voltage 
and  the  simulating  dynamic  response  curve  when  the  load  cur¬ 
rent  changed  from  0A  to  7.5  A,  15  A,  30  A  and  52.5  A  individually. 
The  voltage  curve  went  down  gradually  since  the  current  dragged 
from  the  stack  increased  as  the  load  increased;  at  any  moment, 
more  charges  will  move  and  accumulate  on  the  charge  double 
layer,  which  decrease  the  inner  resistance  of  the  PEMFC  and,  in 
turn,  affects  the  dynamic  response  of  the  fuel  cell.  The  compari¬ 
son  shows  that  the  difference  of  the  average  output  voltage  is  less 
than  0.2  V  when  the  load  increases  to  different  values.  Moreover,  the 
difference  between  the  measurement  and  the  simulation  gradually 
increases  when  the  load  increases,  so  it  can  be  seen  that  the  sim¬ 
ulation  dynamic  curve  totally  agrees  with  the  actual  tested  results, 
and  that  the  model  developed  can  be  used  to  simulate  the  dynamic 
response  of  the  actual  PEMFC.  Furthermore,  it  verifies  the  validity 
of  the  model  developed  and  the  simulation  results. 

To  compare  the  efficiency  of  the  proposed  nonlinear  robust  con¬ 
troller  (NRC),  the  nonlinear  generic  model  control  (GMC)  approach 
[12,13,33]  is  also  implemented  for  the  PEMFC  system.  The  dynamic 
response  of  the  hydrogen  pressure  and  oxygen  pressure  as  the  load 
current  changed  in  step  are  shown  in  Figs.  8  and  9,  where  the  load 
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Fig.  8.  Variation  of  hydrogen  pressure. 
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Fig.  9.  Variation  of  oxygen  pressure. 


current  is  changed  from  0  A  to  92.8  A  for  simulating  the  accelerated 
and  decelerated  process  of  electrical  vehicle.  Fig.  10  displays  the 
absolute  value  of  the  difference  between  the  hydrogen  and  oxygen 
partial  pressures.  It  is  found  in  Fig.  10  that  the  NRC  controller  has 
a  better  transient  response  than  the  GMC  controller.  Generally,  an 
increase  in  the  load  current  causes  a  decrease  in  reactant  pressures 
because  more  fuel  consumption  is  required.  However,  the  flow  rates 
vary  with  the  load  current  in  the  same  way  and  compensate  for  the 
increased  fuel  consumption. 

Figs.  11  and  12  give  the  dynamic  responses  of  the  hydrogen 
and  the  oxygen  flow  rates  under  the  load  current  variations.  It 
is  observed  that  the  oxygen  flow  rate  has  bigger  variations  than 


Fig.  7.  Experimental  results  and  simulation  curves  of  transient  voltage  by  increasing 
the  load. 
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Fig.  10.  Variation  of  pressure  difference  of  hydrogen  and  oxygen. 
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Fig.  11.  Variation  of  hydrogen  flow  rate. 


Fig.  12.  Variation  of  oxygen  flow  rate. 


hydrogen  because  the  oxygen  flow  rate  is  more  sensitive  to  the  load 
current  variation  than  the  hydrogen  flow  rate. 

Figs.  13-15  show  that  the  NRC  controller  is  more  stable  than  the 
GMC  controller  under  the  load  current  variations  for  the  states  Pv, 
an,  Pv, ca.  and  PN2.  Therefore,  the  internal  dynamics  problem  of  the 


Fig.  14.  Variation  of  water  partial  pressure  on  the  anode. 
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Fig.  15.  Variation  of  water  partial  pressure  on  the  cathode. 

nonlinear  robust  controller  is  avoided,  because  applying  the  state 
feedback  exact  linearization  with  the  dynamic  extension  algorithm 
can  make  the  relative  degrees  have  the  same  order  of  the  system  so 
that  the  nonlinear  system  is  global  linearized  and  has  a  satisfactory 
controller,  that  is  not  achieved  in  [12,13]. 

7.  Conclusions 

In  this  paper,  a  MIMO  nonlinear  dynamic  model  of  PEMFC  which 
is  implemented  in  Matlab/Simulink  environment  is  proposed  for 
designing  a  nonlinear  robust  control  by  the  state  feedback  exact 
linearization  with  the  dynamic  extension  algorithm.  By  adding  H0 0 
robust  control  to  the  state  feedback  control  law,  the  steady-state 
error  due  to  parameter  uncertainty  can  be  reduced  and  the  system 
is  guaranteed  to  have  better  robustness.  The  comparison  between 
the  experimental  data  and  simulation  results  shows  that  the  model 
developed  is  validity  and  operable.  The  simulation  results  demon¬ 
strate  that  the  nonlinear  robust  controller  has  better  transient 
responses  to  guarantee  stable  operation  of  the  system  when  the 
vehicle  running  process  is  simulated.  Therefore,  the  proposed  non¬ 
linear  robust  control  strategy  should  be  extended  to  the  design  of  an 
overall  control  scheme  for  PEMFC  system  and  should  be  prolonged 
the  operational  life  of  PEMFC. 
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